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Abstract 

We use the POMWIG Monte Carlo generator to predict the cross-sections for double 
diffractive higgs and di-photon production at the Tevatron and LHC. We find that the higgs 
production cross-section is too small to be observable at Tevatron energies, and even at 
the LHC observation would be difficult. Double diffractive di-photon production, however, 
should be observable within one year of Tevatron Run II. 



1 Introduction 

With the start of Run II at the Tevatron the possibility of detecting a light higgs boson has of 
course been the focus of much attention. One suggested experimentally clean search strategy 
is to look in double diffractive collisions in which the final state contains only intact protons, 
which escape the central detectors, and the decay products of the higgs, in this case two b 
quark jets. The protons would be tagged in roman pot detectors a long way down stream of the 
interaction point, allowing a very precise determination of the higgs mass, using the so-called 
missing mass method [[1|]. With such proton tagging, the available center of mass energy of up 
to 200 GeV means that the preferred higgs mass of around 1 15 GeV is kinematically accessible, 
at least in principle. The missing mass method can be used in the case where higgs production is 
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exclusive, that is the process p + p p + gap + H + gap + p. Unfortunately, recent calculations 
of this rate indicate that it is too low to be observable at the Tevatron [|^]. The inclusive process, 
p + p ^ p + gap + H + X + gap + p, however, was recently estimated by Boonekamp et al. [|3|] 
to be large enough to be detectable. If this turns out to be the case, then the double diffractive 
channel might still provide a relatively clean environment in which to produce the higgs. In this 
paper, we use the POMWIG Monte Carlo to predict the inclusive rate. POMWIG implements 
an Ingelman-Schlein type model of diffraction which has proved successful in describing a 
range of diffractive processes in DIS at HERA [p. We describe the POMWiG physics model in 
Section 2. 

We also consider two other double diffractive processes: dijet production and di-photon 
production. The CDF collaboration has measured double diffractive dijet production rates in 
Run I [[7|], and we can therefore compare our predictions to these data. Ingelman-Schlein type 
models such as POMWiG typically predict rates higher than those observed at the TevatronQ. 
Our strategy in this paper is to attribute this difference to rapidity gap survival probability: 
interactions between spectator partons in the beam protons can destroy a rapidity gap produced 
in a diffractive interaction [|^]. We assume that the factor necessary to scale our dijet predictions 
to data is to a good approximation independent of the hard subprocess, and use it to rescale our 
higgs and di-photon production cross-section predictions. 

We expect the di-photon process to be an interesting process to measure in the first years 
of Run II since, as we shall see, the rates are almost certainly high enough to be observed. 
Di-photon production is dominated by the quark-quark diagram of Figure [T]. This is in contrast 
to the dijet and higgs production processes which are dominated by gluon exchange. The di- 
photon process therefore provides complementary information to the dijet process about the 
structure of the diffractive exchange in a very clean experimental environment. 



POMWiG generates diffractive hard scattering events in the spirit of the Ingelman-Schlein 
model. That is to say it assumes that the production cross-section factorises into a product 
of regge flux factor and corresponding parton distribution function. At present, two regge ex- 
changes are included which we refer to as "pomeron exchange" (iP) and "reggeon exchange" 
(iR). In the current version of POMWiG it is not possible to mix the two exchanges. For partic- 
ular details of the POMWiG default parameterisations (which we use in this paper) see [0]. 

Figure ^ illustrates the mechanism by which we generate higgs boson production in 
POMWiG . In all cases, the higgs is accompanied by "pomeron remnants". This whole sys- 
tem is separated by rapidity gaps from the outgoing hadrons due to the requirement that the 
energy fraction lost by the incoming hadrons, ^, is small (the transverse momentum transfer 
to the outgoing hadrons is also small). More specifically the total cross-section for producing 
higgs bosons in events where the incoming hadrons lose a fraction of their momenta less than 
imax (assuming the top quark is sufficiently heavy) is 



2 Double diffraction in POMWiG 




(1) 



For a recent review see 



and references therein. 
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where r = m\/ s (a/s is the hadron-hadron centre-of-mass energy) and the diffractive gluon 
density for hadron i in pomeron exchange is 

r^max 

9i{x,Q^)= diifip/i{ii)9ip{x/ii,Q'^). (2) 

J X 

fip/ii^i) is referred to as the pomeron flux factor and gjp^PjQ"^) is the gluon density of the 
pomeron, where P is the fractional momentum of the pomeron carried by the gluon. Simi- 
lar expressions exist for the subleading exchange which we generically refer to as "reggeon". 
We use the diffractive parton densities as extracted from diffractive deep inelastic scattering 
measurements by the HI collaboration [|^]. 

The relevant POMWIG subprocesses for di-photon production are shown in Figure [T]. The 
hard subprocess is as in Herwig 6.1 [[T0[]. 




Figure 1 : Di-photon production in POMWiG 




H 



Figure 2: Diffractive higgs production in POMVV^IG 
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Figure 3: The inclusive and double diffractive higgs production cross-sections at Tevatron and 
LHC energies as a function of higgs mass are shown in (a). In (b) the ratio of the double 
diffractive to inclusive higgs production cross-sections is shown. 



3 Results 



In Figure ^ (a) we show the inclusive higgs production and double diffractive higgs production 
cross-sections at Tevatron and LHC energies as a function of nih. We require that ^ < 0.1 
for the diffractive cross-sections. The non-diffractive cross-section was generated using the 
Herwig 6.1 Monte Carlo and the diffractive cross-section was generated using the POMWIG 
default parameters for both the iP and ]R contributions^. The steep fall in the double diffractive 
cross-section with increasing higgs mass at the Tevatron is due to the limited phase space for 
high mass resonance production forced by the requirement that ^ < 0.1. No such dramatic fall 
is present at LHC energies. The effect is visible in Figure |] (b), where we show the ratio of 
double diffractive to inclusive higgs production. The phase space constraints are clearly visible 
in Figures ^ (a) and (b), where we show the ^ and /3 distributions for a higgs mass of 115 GeV. 
P and ^ are forced to high values since the center of mass energy of the two gluons participating 
in the higgs production process s = s^i^2PiP2 > In Figures ^ (a) and (b) we show the 
corresponding distributions at the LHC. Put simply, diffraction as implemented in POMWiG 
looks very much like a hadron-hadron collision at center of mass energy s < ^i^2S, which in 
the case of the Tevatron is of order 200 GeV. It is extremely difficult to produce resonances of 
100 GeV or more given this picture. 

In Table 1 we show the cross-sections for higgs production at the Tevatron for mh = 115 
GeV using two different HI fits for the pomeron structure function, as described in [|^]. Fit 2 is 

^Prior to Herwig 6.3, Herwig overestimated the higgs production cross-section by a factor of two. Actually 
this is quite close to the enhancement typical of NLO QCD corrections [ pi] ] and so we use the uncorrected results 
from Herwig 6.1. 
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process 


a/fb 
(inclusive) 


(j/fb 

r 


(j/fb 
R 


(j/fb 
(IP + IR) 


(j/fb 
(xS') 


H bb, mh = 115 GeV (Fit 2) 


726.5 


0.189 


0.024 


0.21 


0.02 


H bb, nih = 115 GeV (Fit 3) 


726.5 


0.642 


0.024 


0.67 


0.07 


H WW, ruh = 160 GeV 


260 


0.003 


0.0002 


0.0032 




2-f,Et> 12 GeV, < 2 


72477 


265.6 


1009.8 


1275.4 


128 


2-f,Et> 20 GeV, V < 1 


6345 


19.6 


63.2 


82.8 


8.3 


2 7, > 7 GeV, (CDF cuts) 




684 


567 


1251 


125 


2 jets Et>7 GeV (CDF cuts) 




243 ■ 10^ 


87-106 


329.6 • 10*5 


33 ■ 10^ 



Table 1 : Predictions for Tevatron 



process 


a/fb 


a/fb 


o-/fb 


a/fb 




(inclusive) 


r 


R 


(P + IR) 


H bb, rUh = 115 GeV (Fit 2) 


37660 


176 


100 


276 


H bb, TUh = 115 GeV (Fit 3) 


37660 


254 


100 


354 


H WW, ruh = 160 GeV 


19050 


97 


36 


133 



Table 2: Predictions forLHC 



the POMWIG default, whilst fit 3 has a gluon distribution that peaks at high (3. Both distributions 
were found to fit the HI measurement of the diffractive structure function equally well, although 



measurements of diffractive dijet production at HERA suggest that fit 2 is favoured [[12|]. We 



also show the cross-section for nih = 160 GeV, where the dominant decay mode is to WW. 

In order to take gap survival effects into account, we have simulated double diffractive dijet 
production in the kinematic range measured by the CDF Collaboration [^, namely 0.035 < 
< 0.095, 0.01 < < 0.03, —4.2 < 7]jets < 2.4. The jets were found using a cone algorithm 
with radius 0.7. CDF measured the cross-section to be 43.6 ± 4.4 stat ± 21.6 syst nb, a factor 
of approximately 10 lower than our result. We therefore estimate that the gap survival factor 
= 0.1. This compares favourably with theoretical estimates [||,[l3|,[l5]. Our final cross- 
section predictions including our estimate of the gap survival factor are shown in the column 
labelled xS\ 

Table 1 also contains our predictions for double diffractive di-photon production for three 
choices of E^, rf^ and ^ cuts. Due to the dominance of the quark initiated subprocess, the 
diphoton production cross-section is dominated by the poorly constrained subleading exchanges 
at the ^ values which will be available at Run II. We also show the production cross-section in 
the kinematic range used in the CDF dijet measurement where the subleading component is 
significantly reduced. Our prediction of 125 events per fb^^ suggests that this process should 
be observed in the first year of Run II. 

In Table 2 we show our predictions for LHC energies. In this case, we omit the gap survival 
factor since we cannot estimate it from data. We note the recent estimate of ~ 0.01 at LHC 
energies 0]. 
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Figure 4: In (a) the ^ distribution is shown for double diffractive higgs production with higgs 
mass rrih = 115 GeV at the Tevatron. The P distribution is shown in (b). 



4 Summary 



We have used the POMWIG Monte Carlo generator to estimate the double diffractive higgs 
boson and di-photon production rates at the Tevatron and LHC. We conclude that the higgs 
production rate is too small to be observable at the Tevatron. Double diffractive di-photon 
production, an interesting process in its own right, should be observable within the first year 
of Tevatron Run II. The situation at the LHC is less clear, primarily due to the uncertainty in 
the gap survival factors at 14 TeV. If S'^ is not much smaller than that at the Tevatron, then we 
would expect of order 10 events per fb^^ for rrih = 115 GeV, which is a sufficiently high rate to 
see a signal if the background can be controlled. 
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Figure 5: In (a) the ^ distribution is shown for double diffractive higgs production with higgs 
mass rrih = 115 GeV at the LHC. The P distribution is shown in (b). 
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